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Measurement of the production of neighbouring jets in lead-lead collisions 
at = 2.76 TeV with the ATLAS detector 


Abstract 

This Letter presents measurements of correlated production of nearby jets in Pb+Pb collisions at = 

2.76 TeV using the ATLAS detector at the Large Hadron Collider. The measurement was performed using 
0.14 nb“^ of data recorded in 2011. The production of correlated jet pairs was quantified using the rate, 
Rar, of “neighbouring” jets that accompany “test” jets within a given range of angular distance, AR, in 
the pseudorapidity-azimuthal angle plane. The jets were measured in the ATLAS calorimeter and were 
reconstructed using the anti-fc* algorithm with radius parameters d = 0.2, 0.3, and 0.4. Rar was measured 
in different Pb+Pb collision centrality bins, characterized by the total transverse energy measured in the 
forward calorimeters. A centrality dependence of Rar is observed for all three jet radii with Rar found to 
be lower in central collisions than in peripheral collisions. The ratios formed by the Rar values in different 
centrality bins and the values in the 40-80% centrality bin are presented. 


1. Introduction 

Experimental studies of jet production in Pb+Pb 
collisions at the LHC can directly reveal the prop¬ 
erties of the quark-gluon plasma created in the 
collisions. One predicted consequence of quark- 
gluon plasma formation is “jet quenching” that 
refers to the modification of parton showers initi¬ 
ated by hard-scattering processes which take place 
in the quark-gluon plasma [1]. Measurements of 
jet pairs at the LHC provided the first direct evid¬ 
ence of jet quenching [2, 3]. In those measurements, 
the enhancement of transverse momentum imbal¬ 
ance of dijets in central Pb+Pb collisions was ob¬ 
served. Measurements at the LHC of inclusive jet 
suppression [4, 5] and the variation of the suppres¬ 
sion with jet azimuthal angle with respect to the 
elliptic flow plane [6] have shown that the trans¬ 
verse energy of jets is significantly degraded and 
that the energy loss depends on the path length of 
the parton shower in the plasma. These dijet and 
single-jet measurements provide complementary in¬ 
formation about the jet quenching process. The 
single jet measurements are sensitive to the average 
partonic energy-loss while the dijet measurements 
probe differences in the quenching between the two 
parton showers traversing the medium. Those dif¬ 
ferences can arise from the unequal path lengths of 
the showers in the medium or from fluctuations in 
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the energy loss process itself. 

To help disentangle the contributions of these 
factors to the observed dijet asymmetries, the meas¬ 
urement of the correlations between jets that are at 
small relative angles was performed. Neighbouring 
jet pairs include jets originating from the same hard 
interaction, but also jets from different hard inter¬ 
actions. The latter are not of interest in this ana¬ 
lysis, and are subtracted statistically. The remain¬ 
ing neighbouring jet pairs result primarily from 
hard radiation by the parton that occurs early in 
the process of the shower formation. Generally, two 
neighbouring jets originating from the same hard 
scattering should have more similar path lengths 
in the medium compared to the two jets in the 
previous dijet measurement. Therefore measuring 
neighbouring jets could probe differences in their 
quenching that do not result primarily from differ¬ 
ence in path length. More generally, measurements 
of the correlated production of jets in the same par- 
ton shower may provide more detailed insight into 
the modification of the parton shower in the quark- 
gluon plasma beyond the subsequent quenching of 
the resulting jets. 

This Letter presents measurements of the pro¬ 
duction rate of neighbouring jets in Pb+Pb col¬ 
lisions at = 2.76 TeV characterized by the 

quantity Rar introduced in Ref. [7]. The Rar 
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variable quantifies the rate of neighbouring jets that 
accompany “test” jets within a given range of angu¬ 
lar distance, Ai?, in the pseudorapidity-azimuthal 
angle {rj-cj}) plane^, where AR — {Ar\Y + ■ 

Jets were reconstructed with the anti-fct [ 8 ] al¬ 
gorithm using radius parameter values d — 0.2, 0.3, 
and 0.4. In events with test jets with transverse 
energy Et > 70 GeV, further jets are searched for 
within a certain angular distance from the test jet. 

The rate of the neighbouring jets that accompany 
a test jet, Rar, is defined as 


D { TTitest TT'nbr^ 
5-^T / 
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( 1 ) 

where and E^'^ are the transverse energies of 
the test and neighbouring jet, respectively; fVjeT* 
is the number of test jets in a given E^^^ bin and 
is the number of neighbouring jets. Further, 
the Rar quantity was used to define per-test-jet 
normalized spectra of neighbouring jets as 
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Previous measurements of the correlated produc¬ 
tion of neighbouring jets were performed by the D0 
experiment in pp collisions at the Tevatron [7]. The 
measurements by D0 were intended to measure the 
strong coupling constant, Os, and to test its run¬ 
ning over a large range of momentum transfers. The 
measurements presented in this Letter use similar 
techniques and follow notations introduced in that 
measurement. 


2. Experimental setup 

The measurements presented in this Letter were 
performed using the ATLAS inner detector, calor¬ 
imeter, trigger and data acquisition systems [9]. 
The inner detector [10] measures charged particles 
within the interval \ri\ < 2.5. The inner detector 


^ ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The i-axis 
points from the IP to the centre of the LHC ring, and the y 
axis points upward. Cylindrical coordinates (r, (/>) are used 
in the transverse plane, 0 being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of 
the polar angle d as r) = — lntan(6/2). 


is composed of silicon pixel detectors in the inner¬ 
most layers, followed by silicon microstrip detect¬ 
ors and a straw-tube tracker, all immersed in a 2 T 
axial magnetic field provided by a solenoid. The 
calorimeter system consists of a high-granularity li¬ 
quid argon (LAr) electromagnetic (EM) calorimeter 
covering I 77 I < 3.2, a steel/scintillator sampling had¬ 
ronic calorimeter covering It^I < 1.7, a LAr hadronic 
calorimeter covering 1.5 < It^I < 3.2. The hadronic 
calorimeter has three sampling layers longitudinal 
in shower depth and has a A?; x A^i granularity 
of 0.1 X 0.1 for \r]\ < 2.5 and 0.2 x 0.2 for 2.5 < 
|? 7 | < 4.9.^ The EM calorimeters are segmented 
into three shower-depth compartments with an ad¬ 
ditional pre-sampler layer. The forward regions are 
instrumented with copper/LAr and tungsten/LAr 
forward calorimeters (FCal) covering 3.2 < |? 7 | < 
4.9, optimised for electromagnetic and hadronic en¬ 
ergy measurements, respectively. Two minimum- 
bias trigger scintillators (MBTS) counters are loc¬ 
ated on each side at 3.56 m along the beamline from 
the centre of the ATLAS detector. The MBTS de¬ 
tect charged particles in the range 2.1 < \r]\ < 3.9. 
Each MBTS counter is divided into 16 sections, 
each of which provides measurements of both the 
pulse heights and arrival times of energy deposits. 
The zero-degree calorimeters (ZDCs) are located 
symmetrically at z = ±140 m and cover I 77 I > 8.3. 
In Pb±Pb collisions the ZDCs measure primarily 
“spectator” neutrons, which originate from one of 
the incident nuclei and do not interact hadronically 
with nucleons of the other nucleus. 

Minimum-bias Pb±Pb collisions were required 
either to have the transverse energy in the whole 
calorimeter exceeding 50 GeV at the Level-1 trig¬ 
ger or to have a track reconstructed in the inner 
detector in coincidence with ZDC signals on both 
sides. 

Events with high-px jets were selected using a 
combination of a minimum-bias Level-1 trigger and 
High Level Trigger (HLT) jet triggers. The Level- 
1 minimum-bias trigger required a total transverse 
energy measured in the calorimeter to be larger 
than 10 GeV. The HLT jet trigger used the offline 
Pb±Pb jet reconstruction described in Section 4, 
except for the application of the final hadronic en¬ 
ergy scale correction. The HLT jet trigger selected 
events containing a d = 0.2 jet with Ax > 20 GeV. 


^An exception is the third sampling layer that has a seg¬ 
mentation of 0.2 X 0.1 up to \ri\ = 1.4. 
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3. Event selection and data sets 

This analysis used data from Pb+Pb collisions at 
a nucleon-nucleon centre-of-mass energy of = 
2.76 TeV recorded by ATLAS in 2011. It util¬ 
izes data samples corresponding to a total integ¬ 
rated luminosity of 0.14 nb“^. The minimum- 
bias sample was recorded with different prescales 
depending on the instantaneous luminosity in the 
LHC fill. The prescale indicates which fraction of 
events that passed the trigger selection were selec¬ 
ted for recording by the DAQ. The minimum-bias 
trigger recorded an effective luminosity of 7 
Events selected by the minimum-bias trigger and 
the jet triggers were required to have a reconstruc¬ 
ted primary vertex with at least three associated 
tracks each with px > 500 MeV and a differ¬ 
ence between time of pulses from the two sides 
of the MBTS detector of less than 7 ns. A total 
of 51 (14.2) million minimum-bias triggered (jet- 
triggered) events passed the applied event selections 
and were used in the analysis. 

In heavy ion collisions, “centrality” reflects the 
overlap volume of the two colliding nuclei, con¬ 
trolled by the impact parameter of the collisions. 
The centrality of Pb-|-Pb collisions was character¬ 
ized by the total transverse energy meas¬ 

ured in the FCal [11]. The centrality intervals were 
defined according to successive percentiles of the 
distribution ordered from the most central 
(highest to the most peripheral collisions. 

Production of neighbouring jets was measured in 
four centrality bins: 0-10%, 10-20%, 20-40%, and 
40-80%, with the 40-80% bin serving as the ref¬ 
erence. The percentiles were defined after correct¬ 
ing the distribution for a 2 % minimum-bias 

trigger inefficiency that affects the most peripheral 
events, which are not included in this analysis. 

The performance of the ATLAS detector and 
offline analysis in measuring jets in the environ¬ 
ment of Pb-|-Pb collisions was evaluated using a 
large sample of Monte Carlo (MC) events ob¬ 
tained by overlaying simulated PYTHIA [12] hard- 
scattering events onto randomly selected minimum- 
bias Pb-t-Pb events, recorded by ATLAS during 
the same data-taking period as the data used in 
this analysis. PYTHIA version 6.423 with the 
AUET2B tune [13] was used. Three million PY¬ 
THIA events were produced for each of five intervals 
of the transverse momentum of outgoing partons in 
the 2 —2 hard-scattering process, with boundaries 
17,35,70,140,280, and 560 GeV. The detector re¬ 


sponse to the PYTHIA events was simulated using 
GEANT4 [14, 15], and the simulated hits were com¬ 
bined with the data from the minimum-bias Pb-l-Pb 
events before performing the reconstruction. 

4. Jet reconstruction and neighbouring jet 
selection 

Jets were reconstructed within the pseudorapid¬ 
ity interval \r]\ < 2.8. The jet reconstruction tech¬ 
niques described in Ref. [4] were used, and are 
briefly summarized here. The anti-fc* algorithm 
was first run in four-vector recombination mode, 
on A ?7 X A(j) = 0.1 X 0.1 logical towers. The en¬ 
ergies in the towers were obtained by summing en¬ 
ergies of calorimeter cells, calibrated at a scale set 
for electron showers, within the tower boundaries. 
Then, an iterative procedure was used to estim¬ 
ate a calorimeter layer- and 77 -dependent underlying 
event (UE) energy density, while excluding actual 
jets from that estimate. The UE energy was sub¬ 
tracted from each calorimeter cell within the towers 
included in the reconstructed jet. The subtraction 
accounted for a cos 2(j) modulation in the UE energy 
density due to collective flow [ 11 ] of the medium us¬ 
ing a measurement of the amplitude and phase of 
that modulation in the calorimeter. The jet ener¬ 
gies and momenta were calculated via a sum of all 
cells contained within the jets, treating each cell as 
a massless four-vector, using Et values after the 
UE subtraction. A correction was applied to the 
reconstructed jet transverse energies to account for 
jets not excluded or only partially excluded from 
the UE estimate. The magnitude of the correc¬ 
tion was typically a few percent but can be as large 
as 10 % for jets whose energies are fully included 
in the UE estimate. Then, a final 77 - and jet Ex- 
dependent hadronic energy scale calibration factor 
was applied [4]. 

Separate from the calorimeter jets, “track jets” 
were reconstructed by applying the anti-fct al¬ 
gorithm with d = 0.4 to charged particles having 
px > 4 GeV. The track jets were used in conjunc¬ 
tion with electromagnetic clusters to remove the 
contribution of “UE jets” generated by fluctuations 
in the underlying event. The technique is described 
in detail in Ref. [4]. 

In the MG simulation, the kinematics of the refer¬ 
ence PYTHIA generator-level jets (hereafter called 
“truth jets”) were reconstructed from PYTHIA 
final-state particles with the anti-fcj algorithm with 
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radius d = 0.2, 0.3, and 0.4 using the same tech¬ 
niques as applied in pp analyses [16]. PYTHIA 
truth jets were matched to the closest reconstruc¬ 
ted jet of the same d value within AR — 0.2. The 
resulting matched jets were used to evaluate the jet 
energy resolution (JER), the jet energy scale (JES), 
the jet angular resolution, and the jet reconstruc¬ 
tion efficiency. 

The Rar measurement was performed with the 
sample triggered by the jet triggers. The measure¬ 
ment was done differentially in transverse energy 
of the test and neighbouring jets, and in collision 
centrality. Eive intervals 70-80, 80-90, 90- 
110, 110-140, 140-300 GeV and four intervals, 
30-45, 45-60, 60-80, 80-130 GeV were used. Eur- 
thermore configurations where all the Ex bins of 
the test jets or of the neighbouring jets have the 
same upper bound of 300 GeV were also used in 
this analysis. The number of bins and their bound¬ 
aries were chosen to minimize the impact of the lim¬ 
ited number of events in the data while preserving 
the ability to infer the trends in the measured dis¬ 
tributions. Eor each jet radius, neighbouring jets 
are considered if they lie within a specific annu¬ 
lus in AR around the test jet: 0.5 < AR < 1.6, 
0.6 < AR < 1.6, and 0.8 < AR < 1.6 for d — 0.2, 
d — 0.3, and d = 0.4 jets, respectively. The inner 
edge of each annulus was chosen to avoid possible 
overlap of test and neighbouring jets, and the outer 
edge value (< 7r/2) rejects neighbouring jets in the 
hemisphere opposite to the test jet and maximizes 
the number of events. Choosing a maximum Ai? of 
1.6 restricts the pseudorapidity range of test jets to 
jryj < 1.2, yielding approximately 87 x 10^ d = 0.4 
test jets with px > 80 GeV analysed in 0-10% cent¬ 
ral events and 37 x 10^ test jets in 40-80% peri¬ 
pheral events. 

To quantify the centrality dependence of the 
neighbouring jet yields, the ratio = 

7?Afl|cent/7?Ai?|40-80 IS Calculated as the ratio of 
Rar measured in each centrality bin to Rar meas¬ 
ured in the reference (40-80%) bin. 

5. Corrections to neighbouring jet rates 

The raw rates of neighbouring jets include a con¬ 
tribution from neighbouring jets that originate from 
different hard partonic interactions in the same 
Pb-|-Pb collision. This combinatorial background is 
present both in the MG simulation and in the data 
and must be subtracted. It is largest in the low 
bins and it increases with increasing centrality of 


the collision, since the probability for the presence 
of two independent hard scatterings in one Pb-|-Pb 
collision is expected to increase with the number 
of binary collisions. The combinatorial background 
is estimated using the differential yield of inclus¬ 
ive jets (d^iVjet/dpd^dEx) evaluated in minimum- 
bias Pb-|-Pb events. To each event considered a 
weight is assigned such that the event sample ob¬ 
tained from the minimum-bias trigger has the same 
centrality distribution as the sample collected by 
the jet trigger. This estimated background needs 
to be corrected for a geometrical bias present in 
the case where the combinatorial jet overlaps with 
a real neighbouring jet or when two combinatorial 
jets overlap. These biases were removed by apply¬ 
ing a multiplicative correction factor to background 
distributions prior to the subtraction. This multi¬ 
plicative factor was derived from the reconstruction 
efficiency of two neighbouring jets evaluated as a 
function of their angular separation in the annulus. 
In that evaluation, one jet was required to originate 
from PYTHIA’s hard scattering and the other jet 
was required to originate from the minimum-bias 
data in the overlay. The impact of this correction 
on the final subtracted distribution is smaller than 
0.5%. 

The combinatoric jet kinematics may also be af¬ 
fected by the presence of a test jet. To control 
this influence, a study comparing the combinatoric 
jets from the overlay MG events with the same jets 
in the original minimum-bias data was performed. 
This study resulted in an additional correction, in¬ 
dependent of centrality and jet Ex, that decreases 
the combinatorial background by 1.5%. The ±1.2% 
uncertainty on the correction originates from the 
limited number of events and was included in the 
systematic uncertainties. 

In order to account for the effect of the azimuthal 
dependence of jet yields [6], the combinatorial back¬ 
ground was reweighted to take into account the 
measured azimuthal distributions of test jets as well 
as combinatorial jets. The change of the raw sub¬ 
tracted distribution in central collisions and low 
E^’^ bins after the reweighting is at the level of 
1% and decreases with increasing centrality of the 
collision and E^'^. 

The background is subtracted from raw Rar dis¬ 
tributions both in the data and in the MC simu¬ 
lation, allowing an evaluation of the effectiveness 
of the subtraction using the MG simulation. The 
signal-to-background ratio strongly depends on the 
centrality of the collision and E^'^. In 0-10% cent- 
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ral collisions, the signal-to-background ratio can 
be as low as 0.15 for the most extreme case of 
30 < < 45 GeV, and increases to approxim¬ 
ately 0.8 for 60 < < 80 GeV. 

The raw subtracted Rar distributions are af¬ 
fected by the jet energy resolution. The combin¬ 
ation of the jet energy resolution and the steeply 
falling Et spectrum produces a net migration of 
jets from lower E^ to higher ifx values such that 
a jet reconstructed with a given E^'^ corresponds, 
on average, to a lower truth jet Et, . The 

relationship between and E^^ was evalu¬ 

ated in simulated events for the different centrality 
bins and three jet radii used in the analysis. The 
extracted relationships were used to correct for the 
average shift in the test jet energy. No correction 
due to the jet reconstruction efficiency for the test 
jets is needed, since the analysis operates in the 
transverse energy region where the jet reconstruc¬ 
tion is fully efficient. No correction due to jet trig¬ 
ger efficiency is needed either since the plateau of 
the jet trigger efficiency is reached for all test jets, 
except for d = 0.4 jets with E^^* < 90 GeV in the 
0-10% and 10-20% centrality bins. In the region 
70 < i?T < 90 GeV, the jet trigger efficiency is 
above 85%. A systematic uncertainty is applied to 
describe the effect of the lower jet reconstruction 
efficiency. 

The impact of the jet energy resolution, re¬ 
construction efficiency, and angular resolution on 
neighbouring jet yields is corrected for by apply¬ 
ing bin-by-bin unfolding to the raw subtracted Rar 
distributions. For each measured Rar distribution, 
two corresponding MG distributions are evaluated, 
one using truth jets and the other using jets after 
the detector simulation. The ratio of these two MC 
distributions provides a correction factor which is 
then applied to the data. 

The bin-by-bin correction factors are derived 
from the MG simulation where the reconstructed 
jets were matched to the truth jets. To account for 
the impact of the jet angular resolution, the truth 
jet is required to he within a given annulus while 
the reconstructed jet is allowed to fall outside of the 
annulus. 

Examples of jet reconstruction efficiencies for 
neighbouring jets and the bin-by-bin correction 
factors are shown in Fig. 1 for different centrality 
selections and for two choices of jet radii: d = 0.4 
and d = 0.2. Generally, the jet energy resolution in 
central (0-10%) collisions for d = 0.4 jets has com¬ 
parable contributions from UE fluctuations and the 



[GsV] E* [GeV] 

Eigure 1: Reconstruction efficiency (left) of neigh¬ 
bouring jets as a function of the transverse energy 
of the neighbouring jet at Monte Carlo generator 
level bin-by-bin correction factors for 

the distributions of the neighbouring jet rate Rar 
(right) as a function of the reconstructed transverse 
energy of the neighbouring jet E^'^. Plots for two 
different jet radii, d = 0.4 (upper plots) and d = 0.2 
(lower plots) are shown and the transverse energy 
of the test jet E^^^ is required to exceed 90 GeV. 
The four different centrality bins are denoted by 
different markers in each plot. The vertical error 
bars represent statistical uncertainties. 


“intrinsic” resolution of the calorimetric jet meas¬ 
urement. The fluctuations in the UE are approx¬ 
imately two times smaller for d = 0.2 jets than 
they are for d = 0.4 jets. Thus, the distributions 
measured using d = 0.2 jets are far less sensitive 
to the effects of the jet energy resolution, and con¬ 
sequently the resulting bin-by-bin correction factors 
for those distributions exhibit only a modest cent¬ 
rality dependence. The difference in the jet recon¬ 
struction efficiency between the two choices of jet 
radii is also significant - the efficiency for d = 0.2 
jets plateaus around 30 GeV, where it is still rising 
rapidly for d = 0.4 jets. 

The jet angular resolution is determined in MC 
simulation as the standard deviation of the differ¬ 
ence in angles between truth and reconstructed jets. 
In both rj and (f) it reaches 0.008 in 0-10% colli¬ 
sions and 0.005 in the 40-80% centrality bin for 
d = 0.4 jets with E^ = 30 GeV. The angular resol- 
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6. Systematic uncertainties 
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Figure 2: Summary of how the corrections impact 
the dRAn/dE^'^ distribution measured in the data 
in different centrality bins. The dRAR/dE^'^is 
shown for d = 0.4 jets for > 80 GeV and in the 
interval 30 < E!^’^ < 45 GeV. Squares show the raw 
di?Afl/di?y’' prior to the UE subtraction, circles 
show the combinatorial background, triangles show 
the subtracted dRAn/dE^^ prior to unfolding by 
applying the bin-by-bin correction factors, and dia¬ 
monds show the unfolded dRAn/dE^'^. Vertical er¬ 
ror bars on the combinatorial background, raw, and 
subtracted distributions represent statistical uncer¬ 
tainties. Vertical error bars on the unfolded dis¬ 
tribution represent the combined statistical uncer¬ 
tainty from the unfolding and from subtracted dis¬ 
tributions. 


ution improves with increasing jet Ex and reaches 
0.004 (resolution in rj) and 0.002 (resolution in (fi) 
at Ej- = 200 GeV, independently of centrality. 

The impact of these corrections on Rar distri¬ 
butions measured in the data in different central¬ 
ity bins is shown in Fig. 2. The figure shows the 
raw Rar distribution, and the combinatorial back¬ 
ground, subtracted, and final unfolded Rar distri¬ 
butions for d = 0.4 jets with E^^^ > 80 GeV and 
in the lowest E^'^ interval, 30 < E^'^ < 45 GeV. 
The raw Rar distribution increases from peripheral 
to central collisions. The increase of the combin¬ 
atorial background is steeper than the increase of 
the raw distribution. Therefore, a decrease of sub¬ 
tracted Rar with increasing collision centrality is 
observed. The shape of the Rar distribution re¬ 
mains unchanged when the bin-by-bin correction is 
applied. 


Systematic uncertainties in the measurement of 
Rar distributions and their ratios, pr^j^ , arise from 
the uncertainty on the jet energy scale, jet en¬ 
ergy resolution, angular resolution, bin-by-bin un¬ 
folding, centrality, combinatorial background and 
jet trigger efficiency. The impact of uncertainties 
on the jet energy scale, jet energy resolution and 
jet angular resolution was assessed by constructing 
new bin-by-bin correction factors with a systemat¬ 
ically varied relationship between the reconstructed 
and truth jet kinematics. The resulting uncertain¬ 
ties on Rar and were evaluated from their 

changed values obtained with modified jet energy 
scale, jet energy resolution and jet angular resolu¬ 
tion dependencies. 

The systematic uncertainty due to the jet en¬ 
ergy scale is composed of two parts: an ab¬ 
solute, centrality-independent component, and a 
centrality-dependent component. The uncertainty 
on Rar from the jet energy scale uncertainty is 
evaluated by shifting all reconstructed jet trans¬ 
verse energies by ±1 standard deviation of the jet 
energy scale uncertainty. The absolute component 
is determined from the in situ studies of the calori¬ 
meter response; systematic variations of the jet re¬ 
sponse in the MG simulation [16]; and from stud¬ 
ies of the relative energy scale difference between 
the jet reconstruction procedure in heavy-ion colli¬ 
sions, and the procedure used for inclusive jet meas¬ 
urements in 2.76 TeV and 7 TeV pp collisions [17]. 
The magnitude of the uncertainty on the Rar from 
the absolute jet energy scale uncertainty varies from 
2% to 15% as a function of Ex and radius of the 
jet. The centrality-dependent component of jet en¬ 
ergy scale uncertainty [5] was estimated using the 
PYQUEN MC generator [18] to provide a sample of 
jets with modified fragmentation functions consist¬ 
ent with measurements in quenched jets performed 
by ATLAS [19] and GMS [20]. 

The centrality-dependent jet energy scale un¬ 
certainty reaches 1% for 0-10% central collisions 
and less than 0.25% for 40-80% peripheral colli¬ 
sions. The uncertainty on Rar. originating from 
the centrality-dependent component of the jet en¬ 
ergy scale uncertainty increases from less than 1% 
in peripheral collisions to 3% in central collisions. 

The effect of the jet energy resolution uncertainty 
was evaluated by applying modified bin-by-bin cor¬ 
rection factors where the reconstructed jet Ex was 
smeared. The uncertainty on the jet energy resolu- 
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tion is dominated by the uncertainty in the detector 
response. Thus, the procedure used for jet measure¬ 
ments in the 7 TeV pp collisions [16] is used. The 
smearing factor is evaluated using an in situ tech¬ 
nique involving studies of dijet energy balance. The 
systematic uncertainty on Rar due to the jet en¬ 
ergy resolution varies from 1% to 4% depending on 
the jet Et- The centrality-independent jet energy 
scale uncertainty and the uncertainty from jet en¬ 
ergy resolution tend to cancel in the Pr^^^ ratios 
since both the numerator and denominator in the 
ratios are affected to a similar degree by the vari¬ 
ations accounting for the uncertainties. 

The systematic uncertainty on combinatorial 
contributions originates from the previously noted 
uncertainty on the correction factor taking into ac¬ 
count the difference between jets in minimum-bias 
events and combinatorial jets in the overlay. The 
resulting uncertainty reaches ~ 8% in 0-10% cent¬ 
ral collisions at low and rapidly decreases with 
decreasing centrality or increasing 

The systematic uncertainty associated with the 
bin-by-bin unfolding is connected with possible dif¬ 
ferences in the spectral shape between the data and 
the MC simulation. To achieve better correspond¬ 
ence with the data, the simulated jet spectrum was 
reweighted to match the spectral shape in the data 
before deriving the bin-by-bin correction factors as 
described above. Conservatively, the entire change 
in Rar. and pr^^ induced by the use of reweighted 
bin-by-bin correction factors is taken as a system¬ 
atic uncertainty. Typically, this results in 1-2% un¬ 
certainty on Rar- A maximum uncertainty of 5% 
is reached in 0-10% central collisions for Rar eval¬ 
uated for neighbouring jets with Ei^’^ > 30 GeV. 

The uncertainty on the centrality estimation ori¬ 
ginates from the uncertainty on the estimated inef¬ 
ficiency of the minimum-bias trigger. The analysis 
was repeated with modified centrality bins assum¬ 
ing 100% minimum-bias trigger efficiency. The res¬ 
ulting uncertainty is typically smaller than 5% with 
a mild i?T dependence and a negligible centrality 
dependence. 

The uncertainty associated with the jet angular 
resolution is estimated using modified bin-by-bin 
correction factors where the reconstructed jet r] and 
(/) is smeared to reflect a up to ~15% centrality and 
i?T dependent uncertainty on the angular resolu¬ 
tion. The uncertainty on the jet angular resolution 
was estimated by comparing the angular distance 
between track jets and the closest calorimetric jet 
in the data and in the MC simulation. The mag- 



Figure 3: Summary of the relative systematic un¬ 
certainties, in %, on the Rar distributions {6Rar)- 
The systematic uncertainties due to the jet energy 
scale (JES), the jet energy resolution (JER), the 
jet angular resolution, unfolding, jet trigger effi¬ 
ciency, combinatorial contributions and centrality 
are shown for d = 0.4 jets with > 70 GeV in 
0-10% central collisions. 


nitude of uncertainty on Rar from the jet angular 
resolution is smaller than 2%. 

The systematic uncertainty on the jet trigger ef¬ 
ficiency covers a possible bias caused by selecting 
test jets in the region where the jet trigger is not 
fully efficient. This is the case for the d = 0.4 jets 
with E^®®* < 90 GeV reconstructed in the 0-10% 
and 10-20% centrality bins. Eor that Et region, 
the systematic uncertainty was determined as the 
difference between the trigger efficiencies for inclus¬ 
ive jets and jets that were required to have a neigh¬ 
bouring jet. This trigger efficiency difference is less 
than 5% and is independent of the E^'^. 

To avoid statistical fluctuations in the values of 
systematic uncertainties, the weak depend¬ 

ence of the uncertainties is smoothed by a second- 
order polynomial. Systematic uncertainties on Rar 
for d = 0.4 jets are summarized in Table 1 for the 0- 
10% and 40-80% centrality bins. The table shows 
the maximum values of uncertainties for Rar and 
for pr^h- The total systematic uncertainties for 
jets with the other two jet radii are smaller than 
those shown in the table. For the 0-10% centrality 
bin these systematic uncertainties are also plotted 
in Fig. 3 as a function of 
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Figure 4: Rar distributions for d = 0.4 jets (upper) and d = 0.2 jets (lower) evaluated as a function of 
three different columns show Rar distributions evaluated for three different lower bounds on the 
neighbouring-jet transverse energy, > 30, 45, and 60 GeV. The four different centrality bins are denoted 
by different markers in each plot. The shaded bands indicate systematic uncertainties, vertical error bars 
represent statistical uncertainties. The data points and horizontal uncertainties for the 10-20%, 20-40%, 
and 40-80% centrality bins are shifted along the horizontal axis with respect to the 0-10% centrality bin 
for clarity. 



SRar 


0 -10% 

40-80% 

0 -10% 

JES 

15% 

11 % 

7 % 

JER 

4% 

2 % 

2 % 

Angular 

resolution 

2 % 

0.5% 

2 % 

Unfolding 

5% 

2 % 

5% 

Gentrality 

6 % 

6 % 

6 % 

Gombinatoric 

8 % 

<0.5% 

8 % 

Trigger 

5% 

- 

5% 


Table 1: Maximum systematic uncertainties on 
Rar {SRar) and on the ratio of Rar in central col¬ 
lisions and in peripheral (40-80%) collisions pr^j^ 
(^PR&r) foi 'd = 0.4 jets in the 0-10% and 40-80% 
centrality bins. The systematic uncertainty on the 
trigger is applicable only for < 90 GeV. 

7. Results 

The Rar measurement is performed differen¬ 
tially in collision centrality, transverse energy of 


the test jet, E^^^, and transverse energy of the 
neighbouring jet, E^'^, as described Sec. 4. The 
measured distributions are divided into four cent¬ 
rality bins, 0-10%, 10-20%, 20-40%, and 40-80%. 
Fig. 4 shows the fully corrected Rar distributions 
for d = 0.4 and d = 0.2 jets evaluated as a function 
of . The distributions are shown for four cent¬ 
rality bins and three different lower bounds on the 
neighbouring-jet transverse energy, E^’^ > 30, 45, 
and 60 GeV. The shaded error bands on the plots 
indicate the systematic uncertainties discussed in 
Section 6. The Rar distribution exhibits an in¬ 
crease with increasing Sizeable differences 

between the four different centrality bins are ob¬ 
served for all three jet radii. The yield of neigh¬ 
bouring jets is suppressed as the centrality of the 
collision increases. 

To further investigate the centrality dependence 
of neighbouring jet yields, the per-test-jet normal¬ 
ized Et spectra of neighbouring jets defined in Eq. 2 
were evaluated. The resulting differential Et spec¬ 
tra are shown in Fig. 5 for d = 0.4 and d = 0.2 
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0-10% 

II 

O 

2.73 ± 0.23(stat.) ± 0.12{syst.) 

CO 

o 

II 

2.83 ± 0.16(stat.) ± 0.1A{syst.) 

d = 0.2 

2.81 ± 0.15(stat.) ± 0.15{syst.) 


10-20% 

II 

o 

2.85 ± 0.17(stat.) ± 0.13(s2/st.) 

d = 0.3 

2.51 ± 0.15(stat.) ± 0.11(s2/st.) 

d = 0.2 

2.56 ± 0.16{stat.) ± 0.12{syst.) 


20-40% 

II 

o 

2.90 ± 0.12(stat.) ± 0.10(syst.) 

CO 

o 

II 

2.91 ± 0.11(stat.) ± 0.09(syst.) 

d = 0.2 

2.62 ± 0.13(stat.) ± 0.10(syst.) 


40-80% 

II 

o 

3.26 ± 0.15(stat.) ± 0.13(syst.) 

CO 

o 

II 

3.24 ± 0.15(stat.) ± 0.11(s?/st.) 

d = 0.2 

2.99 ± 0.17(stat.) ± 0.11(syst.) 


Table 2: Power-law index n extracted from fits of 
distributions to a power-law function 
cx for > 90 GeV, for four bins in cent¬ 

rality and three jet radii. The resulting fit error 
takes into account combined statistical and system¬ 
atic uncertainties. 


jets and three different lower bounds on the test- 
jet transverse energy, E^^^ > 80, 90, and 110 GeV. 
The same trend of suppression in central collisions 
can be seen as that for Rar evaluated as a function 
of test-jet transverse energy shown in Fig. 4. This 
is a consequence of the steeply falling shape of the 
Et spectra. To better quantify the differences in 
the Et spectra of neighbouring jets, the Et spec¬ 
tra were fitted to a power-law function, oc IjE^'^, 
and the power index was extracted for all three 
choices of jet radius and four centrality bins. The 
results are given in Table 2. The Et spectra meas¬ 
ured in central and peripheral collisions differ in 
the power-law index by approximately two stand¬ 
ard deviations for both the d = 0.4 and d = 0.3 jets, 
suggesting that the Et spectra may be less steep in 
central collisions than in peripheral collisions. 

To quantify the centrality dependence of the sup¬ 
pression of neighbouring jets, the ratios pr^j^ were 
calculated by dividing Rar measured in each cent¬ 
rality bin, except the peripheral bin, by Rar meas¬ 
ured in the peripheral (40-80%) bin. Fig. 6 shows 
PRar evaluated as a function of Ei^^^ and E^'^. 
Some systematic uncertainties cancel in the central- 
to-peripheral ratio as described in Section 6, res¬ 
ulting in distributions that are dominated by 


statistical uncertainties. Ratios are evaluated for 
d = 0.4 jets, which suffer the least from the statist¬ 
ical uncertainties, that are still large. Nevertheless, 
several characteristic features can be observed: the 
PR&r distributions do not exhibit any strong de¬ 
pendence on the suppression factor of 

the most central collisions is at the level of 0.5-0.7 
for all three lower bounds on E^’^] the suppression 
becomes less pronounced with decreasing centrality. 
This is qualitatively consistent with the observation 
of the centrality-dependent suppression of inclusive 
jet yields [4]. In that measurement, the suppression 
of the inclusive jet yields was evaluated in terms of 
the ratio Rcp of the inclusive jet yield in central col¬ 
lisions to the yield in 60-80% peripheral collisions 
spanning the jet pT range of 40-200 GeV. Values 
of i?cp 0.5 were measured in the 0-10% most 
central collisions and exhibited only a weak jet-px 
dependence. 

Contrary to a modest dependence of Pr^j^ on 
the test-jet Et, the pr^,^ evaluated as a function 
of E^’^ suggests a decrease of suppression with 
increasing E^'^. Such a decrease in suppression 
with increasing E^’^ may in fact be expected. The 
jet quenching is generally expected to depend on 
the initial parton energy, but if the splitting hap¬ 
pens such that the two partons have similar en¬ 
ergy, their quenching would likely be comparable 
due to similar in-medium path-length travelled by 
the two partons forming neighbouring jets. Thus, 
in the configuration of E^’^ « E^^^^ the per-test-jet 
normalization effectively removes the impact of the 
suppression. 

8. Conclusions 

This Letter presents results of a measurement 
of the production of neighbouring jets using pairs 
of jets produced at opening angles less than tt/ 2 
in r] — (jj plane. After subtraction of combinator¬ 
ial backgrounds from different hard-scattering pro¬ 
cesses, such jet pairs result from the production of 
multiple jets in the same hard-scattering process. 
As such, it is complementary to previous studies 
of single-jet suppression and dijet asymmetry. By 
probing the relative quenching of a pair of correl¬ 
ated jets in different collision centralities, this meas¬ 
urement opens up the possibility to study the role 
of fluctuations in the jet quenching process. This 
measurement represents a first, exploratory study 
of how the quark-gluon plasma influences the pro¬ 
duction and/or later evolution of the neighbouring 
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Figure 5: The dRAn/dE^’^ distributions for d — 0.4 jets (upper) and d = 0.2 jets (lo-wer) evaluated as a 
function of E^'^. The three different columns show the dRAu/dE^’^ distributions evaluated for three differ¬ 
ent lower bounds on the test-jet transverse energy, E^^^ > 80, 90, and 110 GeV. The four different centrality 
bins are denoted by different markers in each plot. The shaded bands indicate systematic uncertainties, ver¬ 
tical error bars represent statistical uncertainties. The data points and horizontal uncertainties for 10-20%, 
20-40%, and 40-80% centrality bins are shifted along the horizontal axis with respect to 0-10% centrality 
bin for clarity. 


jets from the same parton shower in heavy-ion col¬ 
lisions. 

The jet angular correlations were measured in 
= 2.76 TeV Pb-I-Pb collisions using 0.14 nb“^ 
of data recorded in 2011 by the ATLAS detector at 
the LHC. The measurements were performed using 
jets reconstructed with the anti-fcj algorithm for jet 
radii d = 0.2, d = 0.3, and d — 0.4. The produc¬ 
tion of pairs of correlated jets was quantified us¬ 
ing the rate of neighbouring jets that accompany 
a test jet, Rar, evaluated both as a function of 
test-jet Et and neighbouring-jet E^. A significant 
dependence of Rar on collision centrality is ob¬ 
served in both cases, suggesting a suppression of 
neighbouring jets which increases with increasing 
centrality of the collision. The centrality depend¬ 
ence of the suppression was further quantified us¬ 
ing the central-to-peripheral ratio of Rar distribu¬ 
tions, pr^h . The trends seen in pr^j^ evaluated as a 
function of neighbouring-jet At indicate a decrease 
in suppression with increasing neighbouring-jet At 


which is, however, of limited significance due to the 
limited size of the available data sample. The pr^j^ 
evaluated as a function of test-jet At exhibits a sup¬ 
pression reaching values of 0.5 — 0.7 in 0-10% cent¬ 
ral collisions and does not show any strong depend¬ 
ence on At. This behaviour of the neighbouring 
jet production can be used to constrain the theor¬ 
etical models aiming to describe fluctuations in the 
jet energy loss. 
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